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Abstract 

We present a method devised by Jacobi to derive Lagrangians of any second-order 
differential equation: it consists in finding a Jacobi Last Multiplier. We illustrate the 
easiness and the power of Jacobi's method by applying it to several equations and also a 
class of equations studied by Musielak with his own method [Musielak ZE, Standard and 
non-standard Lagrangians for dissipative dynamical systems with variable coefficients. J. 
Phys. A: Math. Theor. 41 (2008) 055205 (17pp)], and in particular to a Lienard type 
nonlinear oscillator, and a second-order Riccati equation. 

1 Introduction 

It should be well-known that the knowledge of a Jacobi Last Multiplier always yields a La- 
grangian of any second-order ordinary differential equation [10] , [32] . Yet many distinguished 
scientists seem to be unaware of this classical result. In this paper we present again the method 
of the Jacobi Last Multiplier in order to compare the easiness and the power of Jacobi's method 
with that proposed by Museliak et al [H] for the same purpose. We have already presented the 
properties of the Jacobi Last Multiplier in [20]. Some of our papers [IE]-[3T] and the references 
within may give an idea of the many fields of applications yielded by Jacobi Last Multiplier. 

In p3] the authors searched for a Lagrangian of the following second-order ordinary differ- 
ential equation 

x + b(x)x 2 + c(x)x = (1) 

with b(x),c(x) arbitrary functions of the dependent variable x = x(t). After some lengthy 
calculations they found one Lagrangian. In [31] we show that (TjjJ is a subcase of a more 
general class of equations studied by Jacobi [9], i.e. 
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with <p, B arbitrary functions of t and x. We applied Jacobi's method to equation ([I]) and 
showed that many (an infinite number of) Lagrangians can be easily derived. 

In the present paper we show how to obtain many (an infinite number of) Lagrangians for 
the class of equations 

x + f(x)x + g(x) = (3) 

with f{x) and g(x) arbitrary functions of the dependent variable x(t). With the help of the 
Jacobi Last Multiplier standard and nonstandard Lagrangians can be derived without much 
effort. In [15] the author applied his lengthy method to ([3]) in order to obtain at least one 
Lagrangian. 

This paper is organized in the following way. In section 2, we illustrate the Jacobi Last 
Multiplier and its properties [6]- [10], its connection to Lie symmetries [11], [12], and its link to 
the Lagrangian of any second-order differential equations [TO], [32]. We also exemplify Jacobi's 
method with an equation [13J of the class of equations ([1]), i.e.: 

-a + Xx 2 . A . 

x = x— — ^ (4 

and a nonautonomous equation [5] of the more general class of equations ([2]), i.e.: 

■ 2 

x = -— + j. (5) 
x t 

Both examples were not included in [31]. In section 3, we apply Jacobi's method to the class of 
equations ([3]), and show some particular examples such as a Lienard type nonlinear oscillator, 
and a second-order Riccati equation. In section 4, we conclude with some final remarks. 

In this paper we employ ad hoc interactive programs [17J written in REDUCE language to 
calculate the Lie symmetry algebra of the equations we study. 



2 The method by Jacobi 

The method of the Jacobi Last Multiplier [7] -[TO]) provides a means to determine all the 
solutions of the partial differential equation 

n df 
Af = y^ y a i (x 1 ,...,x n )— = (6) 

i=l % 

or its equivalent associated Lagrange's system 

dxi dx2 dx n 



a\ a2 



(7) 
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In fact, if one knows the Jacobi Last Multiplier and all but one of the solutions, then the last 
solution can be obtained by a quadrature. The Jacobi Last Multiplier M is given by 



d(xi,x 2 , ...,x n ) 



MAS, 



(8) 



where 



d(S,ui,V2, ■ ■ ■ ,v n -i) 
d(x 1 ,x 2 , ...,x n ) 



det 



dS_ 

duji 
dx\ 

dx\ 



df_ 

dx n 
du!i 

dx n 



8lo 



n-l 



dx n 



(9) 



and wi, . . . ,u) n —i are n — 1 solutions of (jSJ) or, equivalently, first integrals of Q independent 
of each other. This means that M is a function of the variables (xi, . . . ,x n ) and depends on 
the chosen n — 1 solutions, in the sense that it varies as they vary. The essential properties of 
the Jacobi Last Multiplier are: 

(a) If one selects a different set of n — 1 independent solutions r/i,... ,r/ n _i of equation ©, 
then the corresponding last multiplier N is linked to M by the relationship: 



N = M 



<9(r?i,..., 



Vn-l) 



d(ui,. . .,u} n -i) 
(b) Given a non-singular transformation of variables 

T : (xi, X2, • • • , X n ) ► (^1; %2i ■ ■ 

then the last multiplier M' of A'F = is given by: 

r d{xi,x 2 , ■■■ ,x n ) 



M = M 



where M obviously comes from the n—\ solutions of AF = which correspond to those 
chosen for A'F = through the inverse transformation t _1 . 

(c) One can prove that each multiplier M is a solution of the following linear partial differ- 
ential equation: 

d(M ai 



E 

i=l 



dxi 



0: 



(10) 



viceversa every solution M of this equation is a Jacobi Last Multiplier. 
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(d) If one knows two Jacobi Last Multipliers M\ and M2 of equation (|6|) , then their ratio is a 
solution w of @, or, equivalently, a first integral of ([7]). Naturally the ratio may be quite 
trivial, namely a constant. Viceversa the product of a multiplier M\ times any solution 
uj yields another last multiplier M2 = M\u. 

Since the existence of a solution/first integral is consequent upon the existence of symmetry, 
an alternate formulation in terms of symmetries was provided by Lie [T2]. A clear treatment 
of the formulation in terms of solutions/first integrals and symmetries is given by Bianchi [2j. 
If we know n — 1 symmetries of (jS! / > say 



Ti = }^ij(xx, . . .,x n )d Xj , i = l,n-l, 
Jacobi's last multiplier is given by M = A -1 , provided that where 



(11) 



A = det 



01 
£1,1 



£l,n 



1.1 



l,n 



(12) 



There is an obvious corollary to the results of Jacobi mentioned above. In the case that there 
exists a constant multiplier, the determinant is a first integral. This result is potentially very 
useful in the search for first integrals of systems of ordinary differential equations. In particular, 
if each component of the vector field of the equation of motion is missing the variable associated 
with that component, i.e., dai/dxi = 0, the last multiplier is a constant, and any other Jacobi 
Last Multiplier is a first integral. 

Another property of the Jacobi Last Multiplier is its (almost forgotten) relationship with 
the Lagrangian, L = L(t,x,x), for any second-order equation 



F(t, x, x) 



(13) 



is m 



M 



d 2 L 

dx 2 



where M = M(t, x, x) satisfies the following equation 

d „ , dF 
_ (logM) + _ = . 

Then equation (|13j) becomes the Euler-Lagrangian equation: 

_ d_ fdL\ + d^ = Q 
dt \dx J dx 



(14) 



(15) 



(16) 



4 



The proof is given by taking the derivative of (|16p by x and showing that this yields (|15p . If 
one knows a Jacobi Last Multiplier, then L can be easily obtained by a double integration, 
i.e.: 

L = J (J Md±\ dx + f l {t,x)x + f 2 (t,x), (17) 

where fi and f 2 are functions of t and x which have to satisfy a single partial differential 
equation related to (fT3j) [25] . As it was shown in [25] , f\ , f 2 are related to the gauge function 
G = G(t,x). In fact, we may assume 

dG 

ax 
BG 

h = -ft+Mt'X) ( 18 ) 

where has to satisfy the mentioned partial differential equation and G is obviously arbitrary. 
In [9] Jacobi himself found his "new multiplier" for the class of second-order ordinary differ- 
ential equation^ studied by Euler [3] [Sect. I, Ch. VI, §§915 ff.] ([2]). Indeed Jacobi derived 
that the multiplier of equation ([2]) is given by: 

M = e v{t ' x \ (19) 



as it is obvious from (fT5j) . Consequently in our previous paper [31] we derived a Lagrangian 
of the class of equations ([2]) by means of fH), i.e. 



L = Ig^*)^ + f ^ tj x) + *_ G (t t x) (20) 



with a function of t and x satisfying the following equation: 



^l + e <p(t^) B (t,x) = 0. (21) 
ax 



Equation §3§ is a particular example of the equation considered by Jacobi. In fact from (fT9j) 
and (f2"Uj) we derive: 

m =a^TT' (22) 

and consequently 

L = —4 r - -7^r 7 + ^-G(t, x). (23) 

2(Ax 2 + 1) 2(Xx 2 + 1) dt V ; 17 

This Lagrangian is known [13] . Equation does not possess any Lie point symmetries apart 
translation in t. Therefore Noether's theorem applied to the autonomous Lagrangian L in (|23[) 
yields the following first integral: 

_ ax 2 +x 2 

2(Ax 2 + l)' { > 



lr This is not Jacobi's original notation. 
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Jacobi proved that in the case of a second-order differential equation if one knows a first integral 
and a last multiplier then the equation can be integrated by quadrature (a new Principle of 
Mechanics, indeed) [6], [7]. 

Equation ([5]) is obtained by the symmetry reduction transformation x = u', t = u of the 
third-order equation: 

u'" = -— , (25) 



u 



where u{T) is a function of T. Equation (0) admits an eight-dimensional Lie point symmetry 
algebra and therefore is linearizable. In [25J it was shown that if one knows several (at least 
two) Lie symmetries of the second-order differential equation f|13f) . i.e. 

Tj = Vj(t, x)d t + Gj(t, x)3 x , j = l,r, (26) 

then many Jacobi Last Multipliers could be derived by means of (|12|) . i.e. 

F(t, x, x) 



1 



M„ 





l 


X 




V n 


G n 









dG„ 

dt 
dG m 



x- 



dVn 

dt 
dV m 



(27) 



dt dt 

with (n,m = 1, r), and therefore many Lagrangians can be obtained by means of (|17p . In 
particular, fourteen different Lagrangians can be obtained if the equation admits an eight- 
dimensional Lie point symmetry algebra. We do not look for the fourteen Lagrangians of 
equation ([5]). Instead we use equation (fT9l) to find a Jacobi Last Multiplier and consequently 
a Lagrangian. In fact from (|19p and (|20p we derive: 

x 2 

JIU-y, (28) 



and consequently 



Lag 



x 2 x 2 
2t 



+ -G(t,x) 



(29) 



If one applies Noether's theorem to Lag then the following five first integrals of equation ([5]) 
can be derived: 



Fix 
FI 2 

Fh 
Fh 



2/ 2 
X [X 



2xxy + x 2 t 2 ) 



x 2 x(—x + xt) 



2t 



x 2 x 2 

~w 

x(x — 

XX 

T 



(30) 



xt) 



We underline that the first integrals FI\ and FI4 could not be derived if the gauge function 
G(t, x) was assumed to be equal to zero. 
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3 Equations with space-dependent coefficients 

The class of equations (|3|) has an obvious Jacobi Last Multiplier and therefore Lagrangian if 
the following relationship holds between f(x) and g(x): 

d -a(l-a)f(x) (31) 



dx \f(x) 

where a is any constant 7^ 1. In fact if (|3ip holds then equation (J3j) can be written as 

u + af(x)u = (32) 

i.e. 

/(x) = -il(log«), with ^ = ^ + ^y, (33) 
and thus a Jacobi Last Multiplier for equation ([3]) is 

M = exp ^ /(z)d^ = exp ^~ J d(Iogu)^ = vT 1/a (34) 
and the corresponding Lagrangian is 

L _ + * G(f> B) _ (* + J&L) + - G(t , (35) 

with G(t, x) an arbitrary gauge function. 

We note that this Lagrangian is autonomous and therefore it admits at least the Noether point 
symmetry of translation in t, and consequently the following first integral 

In=(x+ J&L) n/Vl/ ' !lU) (36) 

We would like to remark that because of property (d) of the Jacobi Last Multiplier, then we 
can obtain another Jacobi Last Multiplier M = InM and consequently another Lagrangian 
of equation ([3]). We will not pursue it here any further but one can envision a deluge of La- 
grangians obtained by simply taking any function of the first integral In in (|36H and multiplying 
it for either M in (|34[) or M, and so on ad libitum. 

3.1 Examples 

It is very easy to obtain a Jacobi Last Multiplier and therefore a Lagrangian for the following 
Lienard type nonlinear oscillator: 

x + kxx H x 3 + Xx = 0. (37) 
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In fact we know that M = exp(J* kxdt). Then if one can put the equation in the form 

ill + akxui = 0, (38) 
with a a constant to be determined, i.e. 

Id 

kx = — (logiti) (39) 

a at 

then 

M = u~ 1/a . (40) 

In the case of equation (|3T1) we have 

1 A: 3 

Mi H — kxu\ = 0, with u\ = x H — x 2 + —A, (i.e., a = 1/3). (41) 
3 3fc 



Therefore Mi = u 1 3 and consequent!}!!: 



l '4 + ^ g( ^^ ?tpt|a + ^ g(m) - (42) 

Actually, we can derive another Lagrangian because substituting f{x) = kx,g(x) = ^-x 3 + Ax 
into equation (f3Tj) yields two different a, i.e.: 



^(j) -a(l-a)/ = => 9a 2 -9a + 2 = =^ a M = ^ ( 43 ) 

The case a = | has been considered above. If we substitute a = | into equation (]32p then we 
obtain M2 = u~ 3 / 2 and consequently 



L 2 = V^+ ^G(t, x) = y x + ^x 2 + ^A + A G ( t) x ). (44) 

Equation (|37p admits an eight-dimensional Lie symmetry algebra and therefore is linearizable. 
Moreover one can determine at least twelve more Lagrangians [29]. We note that L\ admits 
one Noether point symmetry while L 2 admits three Noether point symmetries. 
A particular case of ([3T|) is the second-order Riccati equation: 

x + 3xx + x 3 = 0, (45) 

a member of the Riccati-chain [1] . Equation (|45p is linearizable to a third-order linear equation 
by the transformation x = V(t)/V(t), namely (J45]) transforms into V = 0. It also well-known 



We do not take into consideration any nonessential multiplicative constant. 
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that equation (|45p is linearizable by means of a point transformation because it admits an 
eight-dimensional Lie symmetry algebra generated by the following operators: 



ri 


= t 3 (tx 


- 2)d t - t(xt - 


2)(x 2 t 2 + 2 - 2xt)d x 


r 2 


= xt 3 d t 


- (xt-l)(x 2 t 2 


+ 4 - 2xt)d x 


r 3 


= xt 2 d t 


- x{x 2 t 2 + 2- 


2xt)d x 


r 4 


= xtdt - 


- x 2 (xt — l)d x 




r 5 


= xd t - 


x 3 d x 




r 6 


= dt 






r 7 


= td t - 


xd x 




r 8 


= t 2 d t - 


- 2{xt - l)d x 





(46) 



In order to find the linearizing transformation we have to look for a two-dimensional abelian 
intransitive subalgebra [12], and, following Lie's classification of two-dimensional algebras in 
the real plane [12], we have to transform it into the canonical form 

dx i td x 

with t and x the new independent and dependent variables, respectively. We found that one 
such subalgebra is that generated by Ti and Tg = T 2 — Tg. Then, it is easy to derive that 

tx — 1 x 

t 



x(tx - 2) ' 2t(tx - 2) 

and equation ([1311 becomes 



d 2 x 

< 47 > 

We can derive fourteen different Lagrangians by using (|27p and (|17p . Two of these Lagrangians 
admit five Noether symmetries, i.e.: 

and 

Li9 = + i G{t > x) - (49) 

which are derived from 

JLM 56 = - 1 (50) 

and , 

JLM 19 = - ( ^ 2 + t2 ._ 2te + 2)3 , (51) 
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respectively. We remark that JLMqq can be also obtained from ([32]) . In fact equation ([4"5"]1 
can be written as ^ 

u + 3xau = 0, u = x + x 2 , a = —. (52) 

3 

If one applies Noether's theorem to then the following five first integrals of equation (|45p 
can be derived: 

(x 2 t 2 - 2xt + xt 2 + 2) 2 

1 ~ 4(x 2 + x) 2 

(x 2 i 2 - 2xi + it 2 + 2)(x 2 i - x + xt) 

2 ~ 2(x 2 +x) 2 

x 2 + 2x 

73 = 2(x 2 +x) 2 (53) 
(x 2 i — x + xi) 2 

h = 



(x 2 + x 

X 2 

h ■ 



x 2 t — x + xt 



X 2 + X 



while Noether's theorem applied to L19 yields: 
Jni 



x 2 t — x + xt 



2 



x 2 t 2 - 2xt + xt 2 + 2 
(x 3 t - 2x 2 - 2x)xt + (xt 2 + 4)x + (2xt 2 + 3)x 2 
(x 2 t 2 - 2xt + ±t 2 + 2) 2 



(x 2 t — x + xt)(x 2 + x) 

In 3 = ) ,, 2 ivZI , I 2 , 0/2 ( 54 ) 



In 5 



(x 2 t 2 - 


2xt + xt 2 + 2) 2 


1 


> 2 + x) 2 


(x 2 t 2 - 


2xt + xt 2 + 2) 2 




x 2 + 2x 


2(x 2 t 2 - 


- 2xt + xt 2 + 2) 2 



We remark the importance of the gauge function G(t,x). None of the first integrals above, 
apart ^3 and Ins, could be derived if the gauge function was assumed to be equal to zero. 



4 Final remarks 

The Lagrangian which admits the maximum number of Noether point symmetries is that 
obtained by means of the Jacobi Last Multiplier which comes from the determinant (|27p with 
the two solution symmetries, namely the two-dimensional abelian intransitive subalgebra which 
yields the linearizing transformation. We may infer that this is the physical Lagrangian. Also 
the Lagrangian obtained by using the Jacobi Last Multiplier (|34p possesses the maximum 
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number of Noether point symmetries as shown in the case of equation (|3ip . This may explain 
why Lagrangian (|35|) possesses nice physical properties as shown in [3]. 

In the present paper, we do not claim to have been exhaustive in our presentation of the 
application of the Jacobi Last Multiplier for finding Lagrangians of any second-order differential 
equation. Indeed we would like to encourage other authors to apply Jacobi's method to their 
preferred equation. 
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